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ABSTRACT 

We study the expansion of low density cavities produced by Active Galactic Nuclei jets 
in clusters of galaxies. The long term stability of these cavities requires the presence of 
linked magnetic fields. We find solutions describing the self-similar expansion of struc- 
tures containing large-scale electromagnetic fields. Unlike the force- free spheromak-like 
configurations, these solutions have no surface currents and, thus, are less suscepti- 
ble to resistive decay. The cavities are internally confined by external pressure, with 
zero gradient at the surface. If the adiabatic index of the plasma within the cavity is 
T > 4/3, the expansion ultimately leads to the formation of large-scale current sheets. 
The resulting dissipation of the magnetic field can only partially offset the adiabatic 
and radiative losses of radio emitting electrons. 

We demonstrate that if the formation of large-scale current sheets is accompanied 
by explosive reconnection of the magnetic field, the resulting reconnection layer can 
accelerate cosmic rays to ultra high energies. We speculate that the enhanced flux of 
UHECRs towards Centaurus A originates at the cavities due to magnetic reconnection. 

Key words: magnetic reconnection; cosmic rays: ultra-high energies; galaxies: clus- 
ters: general; acceleration of particles 



1 INTRODUCTION 

Active galactic nuclei (AGNs) produce relativistic jets ema- 
nating from the central black hole that shock the intracluster 
medium (ICM), and create over-pressurized lobes. As these 
lobes expand to achieve pressure equilibrium, they form low 
density, high entropy cavities, also referred to as "bubbles" . 
Most bubbles are not expanding supersonically and, thus, 
are in pressure balance with their surroundings. Their rise 
in the cluster potential is controlled by buoyancy, and they 



appear as depressions in the X-ray emissivity (IFabian et al. 



2000 



|Churazov et~aT1|200T) |Birzan et aL]|2004| |Diehl et al 



20081. The jets that create the bubbles contain both plasma 
and magnetic fields. During the active stage the location of 
the jet termination is seen as a radio lobe, while cavities 
are likely to be the outcome of disconnected lobes during 
the quiescent stage. Jets are active for typical timescales 
of ~ 10 7 years, with luminosities of 10 44 erg s _1 . During the 
quiescent stage bubbles rise in the intergalactic medium. 

It is quite likely that the magnetic field plays a cru- 
cial role, however there is no direct measurement yet. Cavi- 
ties survive on timescales much longer than the timescale of 
buoyant rise. Thus, rising bubbles in the ICM must remain 
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coherent, and avoid being shredded by the Rayleigh- Taylor 
instability (RT), Richtmyer-Meshkov instability (RM) and 
Kelvin-Helmholtz instability (KH). On the other hand, sim- 
ple hydrodynamic models "face multiple failures" (Reynolds 



ct al. 2005) when addressing bubble dynamics: in fluid sim- 



ulations the RT and KH instabilities quickly disrupt the 



bubbles on approximately one rise time {e.g. Kaiser et al. 
20051. One of the principal reasons, perhaps, for the fail- 



ure of fluid codes is that the importance of magnetic fields 
has been generally underestimated in the ICM, based on 
the notion that magnetic fields are dynamically subdomi- 
nant. In fact, even dynamically subdommant magnetic field 
can have drastic effects on bubble evolution, stability and 
entropy production. First, magnetic fields in the ICM drape 
around expanding cavities forming a layer of nearly equipar- 
tition magnetic field, which stabilizes the bubble against 



the Kelvin-Helmholtz and Rayleigh- Taylor instability (Dursi 



2007 Dursi fc Pfrommer 20081. Secondly, if the AGN jet 



which blew the bubble carries large-scale magnetic fields, 
it will also stabilize the bubble. Stable magnetic models 



for cavities have been proposed by Benford (20061; Dong 



& Stone (20091; Gruzinov (20101; Braithwaite 



(20101, who 



find the fields inside the bubble and suppress instabilities. 
Examples of magnetic fields in the exterior of the bubble 
that contribute to the stability of the cavities via the mag- 
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netic draping effect have been presented by Lyutikov ( 2006 1 



Nevertheless, it has been argued that unmagnetised bubbles 
can be stable provided a dense shell forms around it during 
the inflation stage ( [Sternberg fe Soker||2008[ ). 

Motivated by the above astrophysical structures we ad- 
dress and solve the following physical problem. Let us as- 
sume that a magnetic structure containing plasma and mag- 
netic field in equilibrium is confined by some external pres- 
sure. This pressure is taken to be constant on the boundary 
of the cavity. As the bubble rises in the gravitational poten- 
tial of a cluster, the external pressure decreases. How will 
the internal structure of the bubble adjust to the changing 
external conditions? We find that the evolution of the cavity 
depends on the adiabatic index of the material: for F = 4/3 
the cavity expands self-similarly keeping the ratio of the ki- 
netic plasma pressure over magnetic pressure constant, while 
for r > 4/3 the structure will be dominated by the magnetic 
field pressure and by kinetic plasma pressure for F < 4/3. 
As a result, for T > 4/3 the expansion of the cavity will 
lead to the spontaneous formation of electric current sheets, 
which are subject to resistive decay of the magnetic field, 
i.e. by forming large-scale reconnection layers. 

To study the evolution of the cavity we use the result 



of previous work by Gourgouliatos et al. (20101 (hereafter 
GBL), who have shown that it is possible to find analyti- 
cal spheromak-type solutions of the structure of a magnetic 
field of spherical geometry, without surface currents, at the 
expense of including some plasma pressure. Those solutions 
are static and stable under MHD instabilities, while the pres- 
sure on their surface is constant allowing them to come to 
equilibrium with a constant pressure environment, without 
any deformation at all. As current sheets occur only when 
there is a discontinuity, those structures have no currents 
sheets and are not susceptible to resistive instabilities, thus 
magnetic flux conservation is a reasonable assumption. 

We study the expansion of the cavity by applying a 
self-similar solution. Substantial work has been done in 
self-similarly expanding magnetic structures (Low 1982 
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(20101 and 



Dalakishvili et al. (2011) have found solutions of force-free 



magnetic fields that expand in a self-similar manner. These 
structures are generalisations of spheromak and Lundquist 
fields, and while having exactly the same geometry, they de- 
pend on time and evolve self-similarly. In this work we apply 
the idea of self-similarity as described in the previous papers 
to the GBL solutions and we study their evolution explic- 
itly in time. For cases where the internal dynamic timescale 
is shorter than the time the external parameters need to 
change substantially we assume quasi-static evolution: a se- 
ries of equilibria solutions which satisfy the changing bound- 
ary conditions. 

In addition to bubble dynamics another astrophysical 
process that is related to the action of AGNs is that of Ul- 
tra High Energy Cosmic Ray (UHECR) acceleration. One 
of the most important results of the newly commissioned 
Pierre Auger Observatory is the discovery of a correlation 
between the arrival direction of UHECRs and AGNs listed 



in the Veron-Cetty & Veron catalog (The Pierre Auger Col 



laboration: J. Abraham et al.||2009| |The Pierre Auger Col- 
laboration et al. 20071. In particular a number of events 



seem to be associated with Centaurus A, the nearst AGN. 
Centaurus A is FR I radio galaxy, with a fa irly low bolo- 
metric radio luminosity of 4 x 10 41 erg s _1 (Alexander & 



Leahy 1987). There is a mildly relativistic jet, detected in 
both radio and at high energies, emitting ~ 5x 10 42 erg s _1 



in ke V - MeV range ( |Kraft et al.||2002[ |Hardcastle et al. 
20031. The orientation of the jet of Centaurus A with re- 



spect to our line of sight is 35° < 9 < 72° (e.g. Skibo et al 



1994). The jet is launched by a central black hole of mass 
Mbh ~ 4 — 10 x 1O 7 M . The galaxy is surrounded by a 
set of AGN-blown cavities, inner (~ 1 — 2 kpc in size), in- 
termediate (~ 10 — 20 kpc in size) and largescale of several 
hundred kiloprsecs, (see, e.g. Israel|[l998[ for review). 



A crucial question is how these cosmic rays are acceler- 



ated. Giannios (20101 has suggested that reconnection can 
accelerate UHECRs up to the essential energies in AGN jets; 



Benford & Protheroe ( 2008 \ have suggested that this process 



takes place in fossil AGN jets through slow reconnection. 
Our suggestion in this paper is that the change of the struc- 
ture of the cavity leads to the formation of current sheets, 
where explosive reconnection is triggered and acceleration 
of UHECRs actually takes place therein. An important is- 
sue with cosmic ray acceleration through reconnection is to 
have slow enough magnetic field dissipation at early stages 
so that a strong field builds, and then a physical mechanism 
to trigger reconnection and accelerate them. The description 
of the cavity evolution we propose provides these properties. 
They start without currents sheets, thus the magnetic field 
gradient is low enough to have a slow dissipation rate, but 
as they expand their plasma /3 drops low enough to create 
a magnetically dominated region and a current sheet. 

The structure of the paper is the following: We start 
from dimensional arguments about the magnetic and plasma 
pressure evolution during the expansion of a magnetized cav- 
ity; then we solve explicitly the partial differential equations 
describing the problem: we find an exact dynamical solution 
and a series of quasi-static solutions which correspond to a 
wider range of parameters. Finally we apply our results for 
the bubble evolution in the acceleration of UHECRs. 



2 DIMENSIONAL ARGUMENTS 

Let us start with a sphere of radius R containing plasma of 
total mass M and and some magnetic flux 9, the plasma 
pressure p is related to the density p through an adiabatic 
index F. The density of the cavity shall be p oc MR~ 3 , the 
magnetic field B oc ^R~ 2 , the magnetic pressure p mag oc 
^ 2 R~ 4 and the plasma pressure p oc p T oc (MR~ 3 ) T . Let 
us consider some change by A_R of the radius of the cavity, 
while the total mass and flux are conserved, and estimate 
the change caused to the plasma pressure, to the magnetic 
pressure and the plasma /3 = p/p mag . 



Ap 
P 



-3r 



AR 
R 

AR 



Pmag R 



(i) 



From the above result we expect that if Y < 4/3 the gas 
pressure will dominate over the magnetic pressure when the 
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sphere expands, whereas if T > 4/3 the magnetic pressure 
will dominate over the gas pressure after expansion, while 
for r = 4/3 the ratio of pressures shall not be affected. 

For r > 4/3, decreasing the plasma /3 implies that at 
some point during expansion, the minimal pressure inside 
the sphere becomes zero. Further expansion leads to the ap- 
pearance of current sheets. 



3 MAGNETIC FIELD STRUCTURE OF THE 
CAVITY 

3.1 The static cavity 

Stability arguments require that the internal structure of 
magnetic fields should be a combination of toroidal and 
poloidal components, which, in order to be in equilib- 
rium with a constant external pressure, must vanish on 
the boundary. These conditions lead to an overdetermined 
mathematical problem: We seek solutions of an elliptical 
partial differential equation, the Grad-Shafranov equation 
I Shafranov 1966 1, with a given value of the function we are 
solving for and its derivative on the boundary. Thus the so- 
lution must satisfy both Newman and Dirichlet boundary 
conditions simultaneously. However, two more functions ap- 
pear in the elliptical partial differential equation: a function 
related to the poloidal current and a function related to the 
plasma pressure. The only constraint they must satisfy is to 
be functions of the flux. As the pressure and the poloidal 
electric current are not determined a priori, they can be 
chosen so that the solution satisfies the over-constraining 
boundary condition. This problem has been solved in GBL, 
and we use this result as the initial state for the rising mag- 
netic cavity. 

In general the magnetic field can be expressed in terms 
of poloidal and toroidal components: 



B = V* x + 2/(*) . 



(2) 



The poloidal flux and current passing through a spherical 
cap of radius r and opening angle 9 are 2-n^{r, 9) and cJ( , 3/) 
respectively. A solution satisfying the Grad-Shafranov equa- 
tion J x B — Vp is 



* = sin 



ci 



a cos(or) 



sin(ar) 



F 
a 



2\ 



(3) 
(4) 
(5) 
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This solution has been studied in detail in GBL, and corre- 
sponds to a linear relation between the physical quantities: 
poloidal magnetic flux, poloidal electric current and plasma 
pressure. Given an appropriate choice of a and Fq the field 
is confined inside a spherical cavity (Figure [T| and the sur- 
face currents can be eliminated, while the whole structure 




Figure 1. Plot of the the poloidal field lines (thin red) and the 
central toroidal field line (thick blue). 



is stable . In the static form of the Grad-Shafranov equa- 
tion the pressure appears through its derivative, while po is 
a constant which is used so that there is no negative energy 
density inside the cavity, thus it has a minimum permitted 
value, but no upper limit. The cavity is a dip in plasma pres- 
sure and density compared to the external medium, whereas 
the magnetic field inside the cavity is higher. The choice of 
po affects /3, the ratio of the plasma pressure over the mag- 
netic pressure in the cavity. The value of j3 is not constant 
within the cavity, we can evaluate an average value 



/ PmagdV 

J pdV 



The value for ft that corresponds to the minimum permitted 
value of po is 0.8, however this is mainly due to the contri- 
bution of the area near the boundaries of the cavity, if we 
stop the integration at 0.9 of the radius then the value of j3 
drops to 0.5, while if we stop at 0.7 the value of j3 is as low 
as 0.2. 

The values for a and Fo are given subject to the condi- 
tions that \&(r c ) = and ^r|r e = 0, where r c is the radius 
of the sphere, these conditions lead to zero fields on the 
surface and a smooth transition from the cavity to the ex- 
ternal medium. For instance, setting r c = 1, and normalising 
the maximum value of VP to unity, we find that a = 5.76, 
ci = -0.134 and F = -24.46. 

It is important to point out that although this solution 
is topologically similar to the spheromak field which is a 
force-free structure inside a sphere, there is a crucial differ- 
ence regarding the surface pressure. The pressure on the sur- 
face of a spheromak is anisotropic. It is zero on the pole and 
reaches a maximum on the equator. Thus a spheromak can- 
not be in equilibrium with a constant pressure environment. 
On the contrary, our solution can be in a stable equilibrium 
with such an environment. 
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3.2 Self-similar expansion 



3.3 Quasi-static evolution 



In this section we allow a — a(t) and we introduce time- 
dependence through self-similar expansion of concentric 
shells with velocity v — —r — r. The evolution of the struc- 
ture is given by the momentum equation, while it must 
satisfy Maxwell equations, the continuity equation and the 
entropy equation. Given that we are working in a non- 
relativistic context we shall use the Maxwell equations keep- 
ing only the appropriate terms. 

(8) 
(9) 
(10) 



V x E 



V • B = 0, 
1 dB 

VxB=- j, 
c 



We shall not consider the dynamical effects of the electric 
field, thus we do not take into account Gauss' law for the 
electric field and in the Ampere-Maxwell equation we shall 
not take into account the displacement current. The electric 
field through the ideal MHD approximation is 



E 



V xB. 



The continuity equation is 
dp 



at 



+ V ■ (pv) = . 



The momentum equation is 



-Vp + 



jxB 



Finally the entropy equation is 
d 



0. 



(11) 



(12) 



(13) 



(14) 



Flux and helicity conservation under self-similar expansion 
lead to the following form for $ 



* = 



n(a(t)7 
a(t)r 



- cos(a(t)r) - F (a(t)r) 2 \ . (15) 



Note that VP = ^(a(t)r,6). The boundary conditions, sim- 
ilarly to the static case, are that the magnetic flux and 
its derivative are both zero at the boundary of the cavity 
r c , while the boundary of the cavity now depends on time 
r c — . The field is expressed in terms of a poloidal and 
a toroidal component. 



B = V* x V(j> + a(i)*V0 , 



(16) 



Constraining our interest to uniform expansion without ac- 
celeration we take from the left-hand-side of the momentum 
equation (131 that a(t) = (vot + ro) -1 , where vo is the ex- 



pansion velocity of the boundary of the cavity. The plasma 
pressure becomes 

F a(t) 4 t> 



P = Po + 



(17) 



From the baryonic mass conservation equation we find that 
the density is 

\3 



p(r,e,t) =p(a(t)r,e)(a(t)r y 



(18) 



where p is a function of th e self similar variable and 6. From 
the entropy equation (141 it is p = Kp F and combining this 

4/3 that po and a arc 



with equation gl8p we find for T 
related through po(i) oc a(t) 4 . 



In the previous section we have found an exact solution for 
the evolution of the cavity. However, this solution constrains 
the value of F and the time dependence of po. In this sec- 
tion we investigate the evolution of the cavity when the ex- 
ternal parameters change slowly compared to the dynamic 
timescale of the cavity itself. Under this assumption we can 
describe the evolution as a series of static solutions which 
are subject to the changing boundary conditions ( |Lynden^| 
Bell 20031, without constraining T, or requiring any time 



dependence for po- Such solutions need to satisfy the right 
hand side of equation ( ] 1 3 p - and as the evolution is slow iner- 
tia terms are negligible; practically we find Grad-Shafranov 
equilibria. 



From equations (171, (18 1 and the adiabatic relation 
(p oc p r ) we find that 



(19) 



Let us express po = p~ou n , where po is constant. Taking 
the logarithm of the above equation we find that there is 
some relation between n and F, which changes as the cavity 
expands and cannot be expressed analytically, except for the 
special case of F = 4/3, where n = 4 and the radius of the 

— 1/4 

cavity is r c oc p which leads us back to the solution of 
section 3.2. 

To find out the approximate form of n when F 7^ 4/3, we 
start from a structure of a given p and then we decrease the 
pressure, demanding that the system remains in equilibrium 
while the mass, the flux and the helicity are conserved. Thus 



the solution shall be of the form of equation (151. Doing so 



for various values of F we find how the radius increases as 
po decreases. In accordance with the dimensional discussion 
of section 2, there are three distinct regimes depending on 
whether F is greater, smaller or equal to 4/3. 

If r = 4/3, P does not change with expansion. This 
is not the case when V > 4/3, since the expansion causes 
a faster decrease of the plasma pressure over the magnetic, 
leading to smaller j3. At some critical point, the gas pressure 
at the minimum becomes zero so that the evolution under 
the self-similar scheme is impossible. At this stage we ex- 
pect that reconnection becomes important. If T < 4/3, as 
the cavity expands its magnetic pressure becomes weaker 
compared to the gas pressure while (3 increases, and fi- 
nally the magnetic field becomes diluted in the background, 
sec Figures |2]l-l(4jl. The inverse process takes place if we 
study the effect of compression instead. From this discus- 
sion we can determine the dependence of r c on po . The gen- 
eral result is plotted in Figure we find r c oc po, where 
n w — (3r) _1 + e n , r appears because of the plasma content 
of the cavity, while e n is related to the contribution of the 
magnetic pressure in the system and its value is such that n 
will tend towards —1/4, with greater contribution when the 
magnetic pressure is important in the cavity. In a similar way 
we evaluate the evolution of P under expansion for various 
values of V. From Figure |6| we find that j3 obeys a power 
law with P oc r ™ , where m — 4 — 3F + e m ; e m is a param- 
eter so that m will tend towards for strongly magnetized 
structures. In general e n and e m will be small compared to 
n and m. Finally we conclude that p oc r c 4 /( 3r ) +1 + £ p w h ere 
e P = (4 — 3r)e n — ^p£m. 
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Figure 2. The pressure for the expanding cavity with T = 1, 
presented as a paradigm of pressure with smaller adiabatic index 
than 4/3. The cavity expands while po drops. As it expands the 
gas pressure dominates over the magnetic one and finally the 
magnetic cavity gets diluted. 
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Figure 3. The pressure for the expanding cavity with V = 4/3. 
The cavity expands while po drops, but it does not become diluted 
nor does it reach zero density at the centre. 

4 ACCELERATION OF UHECR THROUGH 
EXPLOSIVE RECONNECTION 

In this section we shall do some order of magnitude cal- 
culations to show that it is possible to reach the essential 
potentials for the UHECR acceleration. For the numerical 
estimates below we assume that the intrinsic Centaurus A 
jet power is Lj = 10 44 Z/44 ergs s _1 , while during flares it can 
reach 10 45 ergs s _1 . We also assume that the central black 



P 
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Figure 4. The pressure for an expanding cavity with T = 2, 
presented as a paradigm for plasma with V > 4/3. The cavity ex- 
pands while po drops, after some amount of expansion the density 
at the centre becomes zero, while the field there is no longer zero. 
At this stage the field in the centre of the cavity tries to reach a 
force-free equilibrium while the rest of the cavity has a pressure- 
magnetic field equilibrium. This state cannot be described by the 
solution presented in this paper. 
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Figure 5. The dependence of the radius of the cavity r c on po, 
for various values of T, empty circles are for T = 1, asterisks 
r = 6/5, empty squares T = 4/3, cross T = 3/2, solid circle 
r = 5/3 and solid square T = 2. The dependence of r c on po can 
be approximated by a power-law. For T > 4/3 the lines cannot 
be extrapolated to arbitrarily low pressure, after the point that 
corresponds to maximum expansion the minimum pressure in the 
cavity goes to zero and further expansion leads to the formation 
of surface currents. 
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Figure 6. The plasma ft as a function of r c , for various values of 
r, empty circles are for T = 1, asterisks Y = 6/5, empty squares 
r = 4/3, cross r = 3/2, solid circle Y = 5/3 and solid square 
r = 2. The dependence can be approximated by a power law. 
We find that for Y > 4/3 and if j3 ~ 0.8 the minimum plasma 
pressure in the cavity becomes zero and further expansion leads 
to the formation of current sheets. 



hole mass is Mbh ~ 10 s Mq and jet activity lasted for ~ 10 



Myr, close to the estimate of Hardcastle et al. ( 2009 1 of 30 



Myr. We shall show that the weak reconnection throughout 
the evolution of the cavity is not sufficient to reach the avail- 
able potentials. However, as we have shown in the previous 
section, a cavity with initial P > 1 and F > 4/3, after some 
expansion will have parts with zero plasma pressure and 
cannot evolve further in this regime. At this stage some- 
thing must happen, e.g. formation of current sheets, which 
become resistive and form reconnection layers. 



4.1 Unmagnetised cavity dynamics 

The jet ejects material into the intergalactic medium with 
particle density 10 -5 cm -3 . The cavity formed by ejected 
material is confined within a radius r c , further out is the 
forward shock with radius tfs formed from the compressed 
IGM material. From the Sedov- Taylor solution the forward 
shock radius is 

L l/5 t 3/5 



1/5 
PlGM 



60kpc L\^t^ 5 n~_ 



1/5 

5 



(20) 



As the cavity is in equilibrium with the shocked IGM, the 
pressure of the cavity p£ and that of the shocked IGM 
Pigm{jfs /t) 2 are equal. Assuming that most of the energy 
comes in the form of mildly relativistic ejecta L ~ Mc 2 , we 
find for r = 4/3 that 



/ L 7 t 16 \ 
\c 5 pigmJ 



1/30 



30kpc L 4 4 t T 



7/30,8/15 -7/30 



^C^piGM . 

and the ratio of the two radii is 



= 0.5 L 



1/30,-1/15 



TFS 



44 



-1/30 
-5 



(21) 



(22) 



which is a very weak function of the parameters, thus the 
cavity shall be about half the size of the forward shock, 



unless there is an extremely drastic change of the physical 
parameters of the problem. 



4.2 Evolution of a weakly resistive magnetised 
cavity 

Let us assume that a fraction aj ^ 1 of the jet luminosity 
comes out in the form of a large scale magnetic field. Then, 
the jet injects a magnetic flux * at a rate (ajLc) 1 ^ 2 . At the 
same time magnetic flux is destroyed by reconnection at a 
rate f-, where t r is the reconnection timescale. Thus the 
evolution of the flux in the cavity obeys 

■ ^Lcf' 2 ^ 



(23) 



We expect that the magnetic field within the cavity will 
relax to a minimum energy state on an Alfven time scale, 
which as we show below can be very short. Such a minimum 
energy state will include a toroidal-poloidal configuration of 
magnetic field as a purely toroidal or poloidal field is un- 
stable and the co-existence of plasma will make the solution 
described in the previous sections appropriate. As the AGN- 
blown cavity expands, it does work on the shocked IGM, so 
its energy is not conserved; on the other hand, in the absence 
of resistivity the magnetic flux is conserved. Neglecting the 
resistive decay of magnetic flux, the magnetic field in the 
cavity is 

1/30 



B 



^2 



t 2 



7/15 



(24) 



The density in the cavity is very low 

1/10 



Pc = 



(L/c 2 )t 



PlGM 
C 15 t 6 



= 3 x 10"°cm 



-3 7/10 



(25) 



By equating the shocked IGM pressure to p c v^, we can find 
the typical thermal velocity of the cavity particles 



Ut 



— ) 

5 pIGMt 2 



1/20 



= 0.4. 



(26) 



If the jet is ion-dominated, the corresponding ion tempera- 
ture is Ti ~ 4 m p ~ 100 MeV, while for pair-dominated jet, 
the electron temperature is 45 keV. 

When jet ions pass through a reverse shock they are 
heated to Ti ~ m p c 2 ~ 10 9 eV. If electrons are coupled 
to ions effectively, they are heated to the same temper- 
ature. This justifies the use of adiabatic index F = 4/3. 
As the cavity expands, both species cool adiabatically, the 
plasma pressure decreases even faster and the particles be- 
come non-relativistic, leading to V > 4/3, so that the mag- 
netic pressure eventually dominates the cavity. In addition, 
both species can be accelerated non-thermally. 

The Alfven velocity is 



VA 

c 



B 

C^/AlYpc 



( P 7 IGM c35 t 14 



L 7 



1/60 



= 3 



(27) 



Let t]a be the reconnection efficiency, a charge e shall be 
accelerated to energy 



e* = VA V -^Br c = T, A <, i <?'*I*' M t 7 ' 1 V I {& 



n w i r,22 7/20 

2 x 10 eV r/Acrjn_ 5 . 



(28) 



Magnetized cavities and UHECR acceleration 7 



We can neglect resistive decay provided that the reconnec- 
tion timescale, 



1 



L 7 t b 



c35 PlGM 



1/20 



t r ~ r c /(r]AVA) 



is longer than the AGN activity time, t r > t. Thus, 

VI 

t ~ t r > 



(VA V°j) 10/7 C 5 / 2 ^piGM ' 

We should also make sure that the Alfven velocity, which 



(29) 



(30) 



increases with time, Eq. (27 1 does not exceed the speed of 
light 



t < 



e/2 30/14 . 

C ' &j y/PIGM 



(31) 



The most effective acceleration will occur right before the 
resistivity becomes important. Under the constraint that the 
Alfven velocity at this point is smaller than the speed of 
light, Eqns ( |30)31[ ) give <jj = t\a- The resistive time then 
becomes 



and 

VA = 



' L 



Va /7cB/2 VPigm 



L 7 

^agnP'igm 



1/30 



(32) 



10" 2 Llf t- 7/15 nZl /M . (33) 



This value leads to a jet where about 10 -2 of its energy 
budget comes in as magnetic field, thus it is consistent with 
the total luminosity of the system. Using this estimate of 
t)a and the condition <?j — tja, the available potential shall 
accelerate a charge e to energy 



e$ = 



PlGM L AGN 



1/60 



= 2 x 10 is eV. 



(34) 



The main reason the potential ( |34[ ) is fairly small is that 
effective reconnection requires high reconection speeds ($ oc 
v rec ), but if reconnection is too efficient the magnetic field 
is destroyed too early ($ oc B). Still, it might be a viable 
mechanism. First, our order-of-magnitude estimates can eas- 
ily miss a factor of a few; in addition, accelerated particles 
can be of CNO (Centaurus A is located sufficiently close) 
or Fe type. In that case, the acceleration shall take place in 
the area where the field reverses through slow reconnection 
as it has been proposed by |Benford fe Protheroe| ( |2008[ ). 



4.3 Explosive reconnection 



In section |4~2| we have shown that steady reconnection can- 
not lead to sufficient particle acceleration for UHECRs since 
the magnetic field constantly decays, reducing the available 
potential. However, reconnection in astrophysics is rarely a 
steady state process. For instance, solar flares can develop on 
timescales of minutes, while being driven, one might argue, 
on the solar dynamo of 22 years. 

What we have shown in section [53] is that quasi-static 
evolution with the appropriate Y leads to the formation of 
a region with very high magnetization which shall give a 
high value for tja ~ 1 as a strongly magnetized area forms 



inside the cavity. This drastic change leads now to a new 
estimation for the available potential. 



VA T3 

t]A — Br c 
c 



5/4.3/20,7/10 7/20 
VAC L t PiGM 



2xl0 20 L% 20 t 7 7 /10 nZl /2 ° eV. 



(35) 



As the bubble with r > 4/3 expands, at some moment a 
current sheet is formed. We suggest that this serves as a 
reconnection trigger. The ensuing reconnection sets large 



scale flows towards the dissipation region. Equation ( 35 1 



then gives the total electric potential of the reconnection 
layer. This potential is sufficiently high to account for the 
acceleration of UHECRs. 



The Pierre Auger Collaboration et al. ( 2007 1 has found 



that the distribution of UHECRs is not isotropic with a con- 
fidence level of 99%. These observations suggest that there is 
some correlation with AGN sources, especially with Centau- 
rus A, which is a radio galaxy that lies within the GZK hori- 
zon. Centaurus A hosts fossil radio bubbles that are the final 



stages of evolution of the bubbles blown by AGNs (Israel 
1998[ ). In addition, since it is very close, the CRs may com- 
prise of CNO-type nuclei; this decreases the requirements on 
the available potential. 

One may expect two types of acceleration processes in 
reconnection regions. First, there is acceleration by DC elec- 
tric field within the dissipative region itself. This type of 
acceleration is unlikely to produce UHECRs since electrons 
would be accelerated to sufficiently high energies and pro- 
duce pairs, which would annul the electric field. Secondly, a 
Fermi-I type acceleration can occur in the ideal region sur- 
rounding the reconnection layer ( Lazarian fe Vishniac|1999 
Giannios|2010 l. 



5 RADIO LUMINOSITY 

The presence of magnetic fields in cavities leads to syn- 
chrotron radio emission; indeed some X-ray cavities have 



enhanced radio emission (McNamara et al. 2000 Fabian 
|et al.|[2000[ ) which is stronger when the cavity is closer to 
the AGN. During the explosive reconnection stage, it is pos- 
sible to have additional ongoing acceleration within the cav- 
ities, offsetting the effects of adiabatic and radiative cooling. 
However, as we shall show below, this effect averaged over 
a longer time does not change significantly the radio profile 
of the cavities. 

The scenario we consider is the following: When the 
minimum plasma pressure inside the cavity reaches zero the 
ideal expansion cannot proceed any more, current sheets 
form which are susceptible to resistive decay. Resistive de- 
cay proceeds slowly at the beginning, with t\a <C 1, until 
the formation of current sheets with t)a ~ 1. The timescale 
for reconnection at the current sheets can be evaluated from 
equation ( 32 ); it is t r ~ 10 5 . The flux available for reconnec- 



tion is provided with a rate directly related to the expansion 
of the cavity which is longer than the reconnection timescale. 
Indeed, in 10 5 years the cavity shall move with respect to 
the AGN a few kpc or less; assuming a pressure dependence 
in distance from the AGN to be po oc d~ , the expected ex- 
pansion shall be of the order of 2 x 10 -3 of its current radius. 
Demanding that the pressure be non-negative in the cavity 
and repeating the calculations of section 3, we find that the 
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magnetic field must decrease. Assuming mass conservation 
we conclude that about 5 x 10 -4 of the magnetic flux will be 
reconnected. Thus a timescale for reconnection to convert a 
significant fraction of the magnetic field into heat and thus 
radiation is 10 9 years, longer by two orders of magnitude 
compared to the rising time. 

A simple model describing the radio luminosity evo- 
lution is the following. Let us assume that the luminosity 
of the cavity is L ra d = —Nm e c 2 j ra d, where N is the to- 
tal number of emitting electrons in the cavity, and 7 ra( j is 
rate of change of the Lorentz factor, 7, because of radia- 
tive losses. Synchrotron emission depends on the magnetic 
field 'jrad oc B 2 . 7 changes due to three factors: recon- 
nection, which injects energy from the magnetic field in the 
plasma through heating, adiabatic loses to the intracluster 
medium through adiabatic expansion and radiative losses: 
7 = jrec — jexp — "jrad- The timescale of expansion is 10 7 
years, the timescale for radiative losses is 10 s years and as 
we have shown above the reconnection timescale is even 
longer. We conclude that the radio luminosity evolution will 
be mainly affected by the adiabatic expansion, first through 
7 and second through the magnetic field B. These two com- 
bined give the radio luminosity profile shown in Figure Q. 
The radio luminosity will evolve t 3 , where to is 

the rising timescale. The frequency where synchrotron radi- 
ation peaks is io max oc 7 2 B, from the evolution of cavities 
we find that co ma x oc tg 2 , Figure making them harder 
to detect in radio at great distances from the AGN. This 
is consistent with the fact that the cavities that have a sig- 
nificant radio luminosity are the ones close to the AGN or 
directly connected to the jet, which is more evident when 



the evolution is studied in cavities of a given cluster ( Dunn 
et al.|2005 l. 

Thus, although explosive reconnection can produce 
large potentials that accelerate cosmic rays, it does not pro- 
vide an efficient way of converting magnetic energy into heat 
and consequently radio emission. The main reason is that 
while reconnection per se is fast, the flux available to recon- 
nect is limited and determined by the expansion rate which 
is much slower. From a different viewpoint this result can be 
understood in the following terms. If the fraction of the flux 
comprising the reconnection layer over the total flux is 5 and 
the volume of this area is V rec the energy available for radio 
emission is proportional to V^lV^S 2 , whereas the potential 
that accelerates cosmic rays is proportional to V^lV 3 R rec 5. 
In addition even if the volume of the reconnection layer is 
small, one of its dimensions R rec can be large, i.e. of the 
order of the circumference of the cavity, thus the energy 
available for acceleration is 2 x 10 20 eV derived in the previ- 
ous section and therefore sufficient to accelerate UHECRs. 
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Figure 7. The radio luminosity expected by the cavity. While the 
cavity expands by a factor of a few, the luminosity drops almost 
three orders of magnitude, since the evolution of radio luminosity 



is L r 



■ <x t~ 
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to 
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Figure 8. The peak frequency for synchrotron emission. We ex- 
pect that it will move to much shorter frequencies as it moves 
away from the AGN, since its evolution is a; oc t^ 2 



6 DISCUSSION 

In this study we have found analytical solutions for expand- 
ing cavities, both for fully dynamical evolution and also in 
the quasi-static regime. The cavity is confined by some pres- 
sure that decreases with time, but remains constant around 
the cavity. In a more realistic environment and given that 
the cavity size is comparable with the distance the pressure 
of the IGM drops noticeably, we expect a deviation from 



spherical geometry. In addition, the drag force on the top 
surface may cause further deformation. The parametrization 
of such factors can be treated more appropriately through 
a numerical model and is beyond the scope of an analyti- 
cal model such as ours. Another issue is the stability of the 
rising cavity. GBL have shown that static cavities are sta- 
ble using both analytical and numerical methods. In case 
of strong drag forces on the rising cavity, we expect more 
significant perturbations. 
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While the absence of surface currents makes less vulner- 
able to resistive decay, for appropriate values of F current 
sheets will form and their equilibrium will not be described 
by that solution. This allows us to propose that rising mag- 
netized cavities are UHECR accelerators. Expansion leads 
to lower density, and if the gas pressure of the cavity drops 
faster than the magnetic pressure the final outcome will be 
a region containing only magnetic field and confined by cur- 
rent sheets. It is those currents sheets that trigger resistive 
instabilities and allow the available electric potential to ac- 
celerate cosmic rays. Order of magnitude estimations sug- 
gest that this mechanism is viable. 
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